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ABSTRACT: We briefly review the well known connection between classical chaos and 
classical statistical mechanics, and the recently discovered connection between quantum 
chaos and quantum statistical mechanics. 
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Consider a dilute gas of hard spheres in a box with hard walls. Give the spheres some 
arbitrary initial distribution of momenta (and positions). Classically, after a few mean 
free times have passed, we expect that the distribution of momenta will be given by the 
Maxwell-Boltzmann (MB) formula, 

/mb(p) = (27rmA;T)-3/2 exp(-p72mA;T) , (1) 

where the temperature T is given in terms of the conserved total energy U by the 
ideal-gas relation U = ^NkT. 

To see why this should be so, first note that the hamiltonian is simply 

i=l 

where P is a vector with 3N components. Since H takes on the constant value U, the 
allowed values of P form a sphere which we will call the P-sphere. Suppose we now choose 
P "at random." For this to be a meaningful statement, we must have a measure which 
tells us which sets of P's are equally likely a priori. The obvious choice is to assign equal 
a priori probabilities to equal areas on the P-sphere. Then if we choose P at random with 
respect to this measure, the probability that our choice makes an angle between 9 and 
9 + d9 with respect to any particular axis is simply 

f{9) d9 r-. {sin 9)^^-^ d9 

~ {sm9f^-^dcos9 

-(l-COs2^)(3^-3)/2dcOS^. (3) 

If we now identify (2mC/)^/^ cos ^ as, say, the value of pi^ (the z component of the first 
particle's momentum), we find 

fipiz) dpi, ~ (1 - pL/2mt/)(3^-3)/2 dpi, 

~ exp{-pjj2mkT) dpi, , (4) 

where in the second line we have set U = ^NkT and taken the large- A?" limit. Thus we 
have recovered the MB distribution for pi,. Now consider the probability distribution for 
Ply when Pig is fixed; it is given by the first line of (4) with 3N replaced by 3N — 1 (since 
there is one less coordinate when pi, is fixed) and 2mU replaced by 2mU — pi,- In the 
large- A?" limit, we can neglect p\g compared to 2mt/", and we find the MB distribution for 
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Ply. In similar fashion, we get the MB distribution for any n components of P as fong as 
n < TV. 



Now our task is to justify the assumption that equal areas on the P-sphcre are equally 
likely a priori. That such a justification is needed can be seen by considering how we would 
go about filling a real box with a real gas (say, helium). If the box did not already have 
some sort of valve on it, we would install one, and pump the air out through it. Then we 
would close the valve, attach it to a tank of helium with a hose, and open the valve. The 
helium atoms would rush in, moving preferentially in the direction parallel to the hose. 
Thus their initial distribution of momenta would be strongly anisotropic. This is in sharp 
contrast to the prediction of the equal-area measure, which tells us that we will find a 
thermal, isotropic distribution. Clearly, then, the equal-area measure has nothing to do 
with how we put real gases in real boxes, and so we must seek its justification elsewhere. 

That justification comes from Sinai's theorem^, which states that a box of hard spheres 
is a chaotic system. The meaning of this statement in the present context is simple. Start 
off with arbitrary initial momenta and positions; the momenta can be as nonthermal as 
you like. (Actually, we must exclude a set of measure zero; for example, it is possible 
to set up initial conditions such that no two hard spheres ever collide, in which case the 
following discussion obviously does not apply.) Wait a few mean free times, and then note 
the current location of P on the P-sphere. Continue this procedure, keeping track of the 
location of P each time. Chaos implies that this sequence of P's appears to be chosen at 
random with respect to the equal-area measure. 

We are done. Even if we started off with a P representing a strongly anisotropic 
distribution, the next P will appear to be chosen "at random," and so predicts a thermal 
distribution for the individual momenta. 

So much for classical mechanics. What about quantum mechanics? 

Now we have a completely different problem^. The A'^-particle Schrodinger equa- 
tion can always be solved by going to the energy eigenstate basis: H\a) — Uala). The 
hamiltonian is given by (2), supplemented by the boundary condition that the energy 
eigenf unctions 'i/'Q,(X) vanish whenever one of the hard spheres touches a wall of the box, 
or whenever two hard spheres touch each other. The wave function in momentum space 
at time t is then 



where the C^'s specify the initial state. The probability that the first particle has momen- 
tum pi at time t is found by squaring the wave function and integrating over all momenta 




(5) 



a 
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but the first: 



/(Pl,t) = J d'p2...d'pN |V'(P,i)|' 

= ^C:C^e^(^«-^'')*/^ $«/3(pi) . (6) 
In the last hne we have introduced 

*a/3(pi) = J d^P2...d^PN ra{P)MP) , (7) 

which obeys the normalization condition 

J d^pi^af3{pi) = Sal3 ■ (8) 

If we symmetrize or antisymmetrize each V'a(P) on exchange of any two p^'s to reflect 
Bose-Einstein (BE) or Fermi-Dirac (FD) statistics, then /(pi,t) is independent of which 
Pi we choose. 

On physical grounds, we expect that /(pi, t) should be the MB (or BE or FD) distri- 
bution for any time t greater than a few mean free times. It is not obvious how this can 
occur. Consider, for example, what happens if we take the infinite time average of /(pi, t): 



1 \_ 2 

lim - / dt f{pi,t) = Y]\Ca\ ^aa(Pl) 



(9) 



The infinite time average is obviously not something we can actually observe, but theo- 
retically, if anything is going to be thermal, this is it. The problem is that the C^'s are 
essentially arbitrary, so how can we possibly get the MB distribution? 

There is only one way: each $aa(Pi) must individually be equal to the MB (or BE 
or FD) distribution at a temperature Tq, which is given (at least approximately) by the 
ideal- gas relation Ua = ^NkT^- We call this hypothesis eigenstate thermalization. If 
eigenstate thermalization is valid, then (9) will indeed be a thermal distribution as long as 
the uncertainty in the total energy is much less than its expectation value. 

Furthermore, if $a/3(Pi) is always sufficiently small whenever a 7^ then the (3 
terms in (6) will usually make a negligible contribution, and /(pi,t) will be a thermal 
distribution at most times t, without any time averaging at all. However, if the magnitudes 
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and phases of the C^s are carefully chosen, then we can "line up" the $a^(Pi)'s so as to 
get any f{pi,t) that we might want at any one particular time (say, t = 0). Afterward, 
however, as we see in (6), the phases will change in the usual manner; the carefully contrived 
coherence among the various ^*a/3(pi)'s will be destroyed, and we will again find a thermal 
distribution for pi. 

I find this to be a clear and satisfying explanation for the validity of quantum statistical 
mechanics, at least in this particular problem, even without any further evidence in favor of 
it. However, there is more to be said: a very strong case can be made for the two necessary 
ingredients — the thermal nature of $aa(Pi) and the smallness of $a/3(pi) — based on the 
theory of quantum chaos. 

Quantum chaos is the study of quantum systems whose classical counterparts are 
chaotic. The result we will need is known as Berry's conjecture^"^. As its name implies. 
Berry's conjecture is as yet unproved, but there is significant numerical evidence (reviewed 
in [2]) in support of it. 

Berry's conjecture has two parts. Part one says that the energy eigenfunctions of 
a bounded, isolated quantum system which is classically chaotic appear to be gaussian 
random variables, in the sense that 

lim /cZXVa(X + Xi)...^,(X + X,) = Vj(X,^-X,2)...J(X,^_^-X,J. (10) 

pairs 

Here the integration measure is normalized so that / (iX = 1, and the sum is over all 
possible ways to pair up the X^'s; if n is odd the result is zero. Part two says that the 
correlation function J(X) is given by 

J(X) ^ J dP exp{iP-X/n) 6{H(P,X) - U^) , (11) 

where 5{x) is the Dirac delta function, and J(0) = 1. 

It is straightforward to show that Berry's conjecture gives us the two necessary in- 
gredients for quantum statistical mechanics. First, we find that when a ^ 13, $Q,/3(pi) 
is exponentially small in the number of particles A^. More importantly, we find eigen- 
state thermalization: $aa (pi) is given by the MB or BE or FD distribution (depending 
on whether we use nonsymmetric, completely symmetric, or completely antisymmetric 
energy eigenfunctions), plus corrections which depend on the specific energy eigcnfunc- 
tion but which are exponentially small in N. To derive this, the gas must be dilute; 
there are also other, hard-to-compute corrections to ^*aa (Pi) due to the finite radii of the 
hard spheres. We expect these to reproduce the usual hard-sphere corrections to ideal- 
gas behavior, but this remains to be demonstrated. Another important unsettled issue is 
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how high the energy needs to be before (10) is sufficiently accurate. A naive estimate is 
Aq < a, where Aq, = {2nJi'^ /mkTc)^^'^ is a typical particle's de Broglie wavelength, and 
a is the hard-sphere radius. With a in angstroms and m in amu, this condition becomes 
Ta > (SOO/ma^) Kelvin. 

To summarize, the appearance of a thermal distribution of momenta in an isolated, 
bounded quantum system of many particles can only be understood if each energy eigen- 
function individually predicts a thermal probability for the momentum of each constituent 
particle, and if overlaps of different energy eigenfunctions are sufficiently small when one 
particle's momentum is left unintegrated. Both these statements can be derived as conse- 
quences of Berry's conjecture, which is expected to hold only for quantum systems whose 
classical counterparts are chaotic. Thus the well known connection between classical statis- 
tical mechanics and classical chaos is now seen to be mirrored by an analogous connection 
between quantum statistical mechanics and quantum chaos. 

I would like to thank Edouard Brezin, John Cardy, Marty Halpern, and Walter Kohn 
for helpful discussions. This work was supported in part by NSF Grant PHY-91-16964. 

References 

1. Ya. G. Sinai, Ups. Mat. Nauk 25, 137 (1970) [Russ. Math. Surv. 25, 137 (1970)]. 

2. M. Srednicki, |cond-mat / 940303T| , Phys. Rev. E, in press. 

3. M. V. Berry, J. Phys. A 10, 2083 (1977). 

4. M. V. Berry, in Les Houches XXXVI, Chaotic Behavior of Deterministic Systems, 
G. looss, R. H. G. Helleman, and R. Stora, eds. (North- Holland, Amsterdam, 1983). 

5. M. V. Berry, in Les Houches LII, Chaos and Quantum Physics, M.-J. Giannoni, 
A. Voros, and J. Zinn- Justin, eds. (North-Holland, Amsterdam, 1991). 



6 



